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1. Introduction 



More than ten years after its initial discovery [1,2], cosmic acceleration remains an un- 
solved problem. In fact, this phenomenon is so much at odds with conventional particle 
physics and cosmology that a solution might require a complete reformulation of the laws 
of physics governing both very small scales and very large scales. Indeed as we will see, 
the contemporary models trying to explain cosmic acceleration using quantum field theory 
and general relativity fail to provide a convincing framework. In these lecture notes, I will 
not attempt to present the experimental and observational status of cosmic acceleration. 
They have been covered by Jean-Christophe Hamilton in his lectures. I will only try to 
provide some indications about cosmic acceleration. The number of important topics has 
grown enormously in the last ten years, so much that I will only be able to cover a limited 
number of them. 

Cosmic acceleration was first observed using type la supernovae and their Hubble 
diagram, i.e. the redshift vs luminosity distance. The result is purely kinematical and 
stipulate that the acceleration parameter of type la supernovae q = — ^§ is negative in 
the recent past of the universe. This implies that distances measured according to the 
Friedmann-Robertson- Walker metric (FRW) are increasing fast a > 0: 



ds 2 = -dt 2 + a 2 (t) 



,h ' 2 . r 2 {d9 2 ■ sin 2 0# 2 ) 



1 — kr 2 



(1.1) 



Notice that a has a dimension here and the parameter k = 0, ±1 is the reduced curvature. 
Spatial sections are open, closed or flat depending on k = — 1, 1 or k = 0. Another period 
of accelerated expansion seems to have also existed in the early universe and bears the 
name of cosmic inflation [3]. Here the late time acceleration started when z ~ 1 where the 
redshift is defined by 

l + z=^ (1.2) 
a 

and a subscript denotes the present value of a cosmological quantity. 

Understanding the observation of cosmic acceleration requires a theoretical framework. 
In the last century, cosmology [4] has been very successful in describing the evolution of 
the universe using two fundamental assumptions. The first one is: 

The Universe can be described using the general theory of relativity. 

General relativity [5] has been tested in the solar system and beyond, notably in extreme 
astrophysical situations such as binary pulsars. So far, there is no reason to doubt its 
validity up to cosmological scales. Another principle is usually necessary to simplify the 
analysis of the universe as a whole (the cosmological principle) : 

The universe is homogeneous and isotropic on large scales. 

Of course, the universe is lumpy on small scales. Nevertheless, the appearance of small 
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scale structures can be understood as resulting from the growth of initial inhomogeneities. 
The cosmological principle is also a very robust hypothesis. 

General relativity relates the energy content of the universe to its geometry. Obser- 
vations have revealed four different types of energy in the universe. Ordinary matter is 
described by the standard model of particle physics (both baryons and leptons) and is 
responsible for the existence of stars and galaxies. Radiation in the form of photons is the 
best probe we have to observe cosmic phenomena. Neutrinos are elusive particles which 
participate in the radioactive phenomena leading to the creation of the elements (Big Bang 
Nucleosynthesis). Finally, a host of phenomena including the rotation curves of galaxies 
seem to lead to the existence of exotic particles in the form of dark matter. These four 
types of energy are enough to describe most of the history of the universe. Unfortunately, 
they cannot account for a period of cosmic acceleration. The dynamics of the scale factor 
are governed by 

^ = -^(p + 3p) (1.3) 

Ob O 

where Gn is Newton's constant, p = Ylt=i Pi the total energy density of the four types of 
energy and p = Y2i=\Pi the total pressure. Each fluid has an equation of state Wi such 
that pi = Wipi. Ordinary matter and cold dark matter have a vanishing equation of state 
Wi = while radiation and neutrinos have iWj = 1/3. Of course, this implies that a < 0. 

A negative pressure is not enough to guarantee acceleration. One must impose that 
the total equation of state p = wp satisfies w < —1/3. It happens that a fluid with this 
property was introduced by Einstein in order to guarantee the existence of a static and 
spherical universe comprising only ordinary matter. Indeed the dynamics of the universe 
can be equivalently described using the Friedmann equation 

ff2 = 8^_* (14) 

where H = - is the Hubble rate. Imagine that there exists another fluid with p = —p = 
8n Q N on top of ordinary matter. Then the k = 1 Friedmann equation yields: 



2 8ttG n A 1 
H =— pm+ 3-^ (L5) 



where p m is the matter density in the universe, together with 

a At:G n A 

- = 5— Pm + (1-6) 

a 6 3 

A static solution exists with a = Unfortunately, the discovery by Hubble that the 
universe is expanding rules this model out. 

A lesson can be learnt from Einstein's static universe. Indeed, the existence of a 
cosmological constant is a key ingredient in order to obtain a repulsive behaviour in general 
relativity. For Einstein's universe, the negative pressure due to the cosmological constant 
counterbalances the gravitational attraction. If this is not the case, then the repulsive 
nature of the cosmological constant will lead to the acceleration of the universe. This is 
the simplest explanation of cosmic acceleration. 
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The acceleration of the universe can be easily formulated in terms of the energy content 
of the universe. Acceleration requires that the cosmological constant must play a dominant 
role. Indeed this conclusion follows from the type la measurements and the WMAP results 
on the Cosmic Microwave Background (CMB) [6]. Let us denote by 



where the critical energy density is 



fii = - (1.7) 



(LB, 



87tGjv 



We have defined the energy density 1 pa = 8w q , and we consider that there are now 
five fluids including the one corresponding to the cosmological constant. The Friedmann 
equation can be rewritten 

where f2 = Yli=i Observations of the CMB tells us that the universe is almost spatially 
flat fi = 1. This follows from the location of the acoustic peaks in the CMB spectrum. Using 
general relativity and the cosmological principle, one can relate the acceleration parameter 
to the fraction of cosmological constant and matter £l m (we neglect radiation and the 
neutrinos) 

q = ^-n A (1.10) 



With these two relations one deduces that Qa does not vanish and even dominates now 



ft m ~0.3 O\~0.7 (1.11) 

The universe is accelerating now, and the energy content of the universe is dominated by 
a pure cosmological constant. This result which was first deduced about 1998 has been 
greatly refined since then using CMB data, large scale structures and type la supernovae. 
Part of these lectures will be devoted to understanding the physics behind the cosmological 
constant. 

Observations give that 

p A « 10" 48 (GeV) 4 (1.12) 

This density is approximately 10~ 29 g/cm 3 . It is an extremely small scale which is only 
matched by the neutrino mass m v ~ 1CT 3 eV. It is extremely difficult to justify the existence 
of such a small cosmological constant using quantum field theory alone [7]: this is the 
cosmological constant problem. Postulating the existence of a small cosmological 
constant is so problematic that other possibilities have been considered. They fall within 
two categories. The first one amounts to introducing a new type of matter called dark 

x We will freely use the equivalence 8tyGn = k\ = rripf where rapi « 2.10 18 GeV is the reduced Planck 
mass. 
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energy whose role mimics a pure cosmological constant. The Einstein equations governing 
the evolution of the universe 

Rpu - \g^R = 8ttG n T^ (1.13) 

where R^ u is the Ricci tensor, R the Ricci scalar and T^ u the energy- momentum tensor of 
the four types of energy, is supplemented with a new energy-momentum tensor due to the 
dark energy component of the universe 

1 



- ~g»vR = 8ttGn {T^u + T° E ) (1.14) 



Several examples of dark energy will be presented. Another possibility concerns gravity 
itself. Indeed, one could envisage that gravity which has been tested up to galactic scales 
must be modified on larger scales. This could be formulated with a modified Einstein 
equation 

/ (R, R^iv-i Ra.fip(j)iiv = 

8nG N T^ (1.15) 

where is a tensor involving all the components of the Ricci tensor R a j3 pa which reduces 
to the Einstein tensor = R^ — \g pv R on small scales (up to the galactic scales). As 
we will see, this approach is fraught with difficulties. 

These models are easier to describe using the Lagrangian formulation of general rela- 
tivity 

The first term is the Einstein-Hilbert action. The second term involves the cosmological 
constant. Notice that it is a simple additive correction to the Einstein-Hilbert action. In 
fact, this additive constant plays the same role as a vacuum energy. This degeneracy is the 
origin of the cosmological constant problem. 

Using the action principle, it is easy to understand how to incorporate a dark energy 
component, one considers the action 

S = J d 4 x V^(R - £de) (1.17) 

where £de is the dark energy Lagrangian which reduces to a constant on cosmological 
scales. We will give many examples of dark energy Lagrangians. Modifying gravity can be 
easily implemented too, it only amounts to altering the Einstein-Hilbert term 

S = y^q^ J d 4 x ^g~h(R, R^R^pa) (1.18) 

where h is a scalar function. Different forms of h will be presented too. 

So far, we have assumed that the Einstein equations need to be modified in order to ac- 
commodate the acceleration of the universe. In fact, there is another appealing explanation 
which only requires to assume that the cosmological principle is violated. In particular, 
if we lived in a large void, the acceleration of the universe would be only apparent and 
due to matter depletion in the void. Of course, this would imply that observers have a 
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special place in the universe, contrary to the Copernican principle. In the course of these 
lectures, we will encounter other even more drastic routes which have been followed in 
order to understand cosmic acceleration; for instance we will discuss DGP gravity [8] and 
unimodular gravity [9]. 

The acceleration of the universe is a rich subject and I cannot cover all the diverse 
solutions which have been proposed. The following choice of topics is personal. More 
thorough coverage of the subject can be found in excellent review articles [10-12]. Citations 
will unfortunately be parsimonious, mainly general references where more information can 
be found. 

In a first part I shall recall the status of the cosmological constant problem. I shall 
then describe dark energy and its siblings the scalar-tensor theories including f(R) gravity. 
Another chapter is devoted to modified gravity in the infrared. Finally, the possibility of 
a violation of the Copernican principle will be analysed briefly. 

2. The Cosmological Constant 

Cosmic acceleration may be explained by the presence of a cosmological constant supple- 
menting Einstein's general relativity. From the observational point of view, this is by far 
the simplest and most economic explanation. From a fundamental point of view, the cos- 
mological constant term is puzzling. As mentioned, a cosmological constant plays the same 
role as a vacuum energy. Quantum theory has taught us that the vacuum is not empty 
space: it is full of vacuum fluctuations in the guise of particle-antiparticle creations over 
extremely short time scales. These fluctuations have had a nice experimental confirmation 
in the Casimir effect whereby two metallic plates attract each other under the influence of 
electromagnetic quantum effects. As a result, one should think of the cosmological constant 
as an energy density comprising two terms 

p A = A m| 1 + y (2.1) 

where Vq contains all the quantum fluctuations due to the particle physics vacuum and the 
term Aqtu^ is called the bare cosmological constant. The net result given by observations 
is that pa is very small. How can this be? 

2.1 The cosmological constant problem 

Vacuum contributions to the cosmological constant arise from the creation and annihilation 
of particle-antiparticle pairs. Each particle i of mass rrii contributes 

Vb = X>VS, m = -^m\ (2.2) 

i 

where Cj = 0(1). Taking into account known particles such as the gauge bosons Z and 
W, one gets a contribution of order M| which is sixty orders of magnitude larger than 
the observed p\. Facing such a major difficulty, two logical possibilities can be envisaged. 
First of all, there could be a cancelation effect between Vq and Aom P1 . This would require 
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a very precise fine tuning of sixty orders of magnitude. Another possibility is that both Vq 
and Aompj are small. 

Large cancelations are highly unnatural in quantum field theory. For this reason, 
an almost exact cancelation between Vq and Aom pl must have another origin. A clue 
is given by Weinberg's bound on p\. Suppose that observations did not indicate that 
p\ ps 10 _48 (GeV) 4 . Could we infer a reasonable value for p^J This is the question which 
was positively answered by Weinberg [7]. Indeed as soon as a cosmological constant term 
starts dominating the dynamics of the universe, structures stop growing. This implies that 
galaxies must have formed before the start of an accelerated period. This leads to a bound 

PA < 500p m0 (2-3) 

where p m o is the present matter density. In fact this entails that p\ < 10 _46 GeV 4 , i.e. 
a very stringent bound. Thus a strong bound on the cosmological constant results from 
the existence of galaxies only. This reasoning can be extended and stated as an anthropic 
result: the existence of observers is only compatible with a small cosmological constant. In 
a sense, the mere fact that we observe the universe imposes that the cosmological constant 
must be small. As we will see in the next section, this interpretation can be put on firmer 
ground when a large number of possible vacua exist, and if the probability of existence of 
these vacua is evenly distributed, one could argue that there is a non-negligible probability 
that our universe is atypical and happens to accommodate a large cancelation between 
quantum effects and the bare cosmological constant. 

Within the realm of quantum field theory, the previous cancelation is unnatural. To 
this end, it is useful to recall a few fact about effective field theories [13,14]. Let us consider 
a quantum field theory describing high energy phenomena up to a scale E. In this high 
energy theory, a bare cosmological constant is present Pa{E). Let us now consider the 
physics at lower scales p < E. When lower energies are probed, particles of intermediate 
masses p < rrii < E cannot be produced in experiments. They can only appear in the 
vacuum fluctuations. Hence at the lower energy scale p, the physics can be described 
by an effective field theory comprising only particles whose masses are mi < p and a 
cosmological constant term 

PA (p)= PA (E)+ ^2 m t ( 2 - 4 ) 

fi<rrii<E 

In this context, the bare cosmological constant Ao = n\p\{E) is a parameter depending on 
high energy physics whereas the quantum fluctuations involve particles which have been 
integrated out. Describing the acceleration of the universe requires to consider low energies 
well below the electron mass p < m e . Viewed from this angle, an exact cancelation between 
both terms is hard to justify and would require very special properties. It is more natural 
to assume that both pa(E) and the quantum corrections are small. We will see that this 
is also extremely difficult to realise. 

2.2 The Landscape 

Recently, anthropic ideas and string theory [15] have led to the landscape picture whereby 
a large number of vacua could explain the existence of a small cosmological constant. 
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Although the string picture is too complex to present here, we can extract some of its key 
ingredients in a 4d context and a conventional field theory. Let us consider the action 



S = f d 4 xv^(^ - PA K) -^X> 2 ) ( 2 - 5 ) 



where we have taken into account p A (m e ) at the electron mass for the sake of the argument. 
The action involves a collection of J four forms F pi/pa . The equations of motion for the 
four forms yield F pvpu = riiqe^ vp(T where rii is an integer (this is analogous to the Dirac 
charge quantisation). This immediately leads to an effective vacuum energy 

2 J 

PA = pA(m e ) - — n\ (2.6) 
i=i 

Each choice of the integers ni correspond to a different vacuum. In fact, if we consider the 
state vector {n\ . . . nj} as representing one type of universe with its particular cosmological 
constant, can we arrange the sum ^- Yl%=i n i *° cancel PA{rn e )l This would select the type 
of universe we live in. 

The number N of vacua in the interval [n 2 , n 2 + dn 2 ] with n 2 = J2i=i n i ^ s related to 
the area of a J-sphere of radius n, i.e. 

dN _ (2ny/ 2 r 2 



dn 2 2r(J/2) 



n J - 2 (2.7) 



Increasing iV by one unit requires 5n 2 = this is the typical spacing between values of 
n 2 in the space of all vacua. This corresponds to a step in the cosmological constant 

q 2 2T( J/2) 2 _j 
SPA = -J^j^ n (2 - 8) 

Typically, one would like to cancel />A(m e ) which requires to cancel a large number 

n> = 2 f (2.9) 

This cancellation cannot be exact as the right-hand side is not always an integer. This 
is almost exact if the interval between each step is very small, what is left over being 
interpreted as the energy density driving cosmic acceleration. The step size is simply 

8 PA _2V { J,2) n _ J 



PA (m e ) {2it)J/ 2 

This ratio measures the amount of fine tuning in order to cancel the cosmological constant 
approximately and gives a bound on the leftover cosmological constant which drives the 
cosmic acceleration. It decreases with an increasing number of forms J. Canceling sixty 
orders of magnitude requires ) = 10~ 60 which can easily be achieved with a large 
number of forms. 
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This toy model is a good illustration of the landscape paradigm. Indeed a large number 
of four forms can almost cancel the cosmological constant, leaving a tiny remainder whose 
presence leads to the current cosmic acceleration. Although the adjustment is fine tuned, 
there is a non- vanishing number of vacua specified by the integers rij which can satisfy these 
constraints. If all the possible vacua can be populated, the anthropic principle specifies 
that we simply happen to live in one of these particular universes. 

2.3 Weinberg's theorem 

So far we have not tried to have a dynamical understanding of cosmic acceleration. In- 
deed the landscape picture just postulates that the overall effect of the bare cosmological 
constant and quantum fluctuations can be almost completely canceled by negative contri- 
butions due to many four forms. From the point of an effective field theory, this is not 
satisfying and appears as a last resort explanation. A more natural explanation would 
be the existence of a symmetry principle which would guarantee that the cosmological 
constant vanishes. This symmetry may be slightly broken resulting in a tiny cosmological 
constant whose value would drive the cosmic acceleration. This is the technical meaning of 
natural in field theory. A theory is natural if by taking a parameter to zero in a Lagrangian, 
the degree of symmetry of the theory increases. We will present two such symmetries and 
analyse their drawbacks: scale invariance [7, 13, 14] and supersymmetry. 

The cosmological constant pa has mass dimension four. This implies that a change 
of scale x — > Xx rescales p\ — > A _4 pA- Of course, if the field theory describing cosmic 
acceleration is scale invariant, the resulting cosmological constant must vanish 

PA = (2.11) 

The simplest field theories which could model out cosmic acceleration are scalar field the- 
ories. Indeed scalar fields <pi do not carry Lorentz indices and can then acquire a vacuum 
expectation value (vev) without breaking Lorentz invariance. Scalar field theories can be 
specified by their interaction potential V((j)i). In particular we always include the bare cos- 
mological constant as a constant term in the potential. Vacua of these theories minimise 
the potential 

Wfc=<fc> = (2.12) 

We assume that the potential V(4>i) is the low energy effective potential valid for scales well 
below the electron mass. Hence it captures the effect of integrating out massive particles. 
The vev of the effective potential is then interpreted as the effective cosmological constant 

PA = V(< & >) (2.13) 

Now let us assume that the low energy effective action is classically scale invariant. This 
implies that the potential can be expanded as 

v(4>i) = Xijki(f>i<Pj4>k(t>i ( 2 - 14 ) 

in terms of monomials of dimension 4. Each scalar field has dimension one. We can always 
single out one field <j>i and write 

Vfa) = <t>\V{ Zi ) (2.15) 
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where = (pi/ (pi and V(zi) = V(l, Z2 • • • zjv)- Vacua of this theory satisfy 

3^ = 0,^ = 0, (2.16) 

i.e. they have a vanishing cosmological constant. The potential V is a function of (N — 1) 
variables. The existence of a minimum with a vanishing cosmological constant is only 
possible if there exists a relation amongst the coupling constants 

g(\ ljk i) = (2.17) 

Notice too that the vev of fa is not specified, it is a Goldstone direction corresponding to 
the spontaneous breaking of scale invariance by the non-trivial vacuum < Z{ >. The origin 
where all the fields fa vanish is also a vacuum of the theory. In fact, writing < fa >= fa < 
Zi >, one can see that there is a flat direction of vacua in field space parameterised by fa. 
Along this flat direction, the cosmological constant vanishes identically. 

The scalar field fa can be used to give masses to fermions thanks to Yukawa couplings 

£yuk = Vijkfaipjipk (2.18) 

leading to fermion masses rrijk = yijk < (pi >■ Apart from the origin, all the vacua along 
the flat direction lead to non- vanishing masses. Although the effective action at low energy 
is obtained by integrating out massive particles and only contains massless fields before 
spontaneous breaking of scale invariance, quantum corrections are still present and lead to 
a logarithmic running of the couplings constant Xijki- This scale dependence implies that 
generically 

9{\jM)*U (2.19) 
unless another symmetry is present and guarantees that the relation between the coupling 
constants is invariant under the renormalisation group evolution. If this is not the case, 
the vacuum < zi > is lifted by quantum corrections. Of course, there always remains the 
trivial vacuum where all the fields vanish. Unfortunately this vacuum cannot give masses 
to fermions and is therefore not suitable to describe our universe. 

Weinberg's theorem is very powerful as it rules out most attempts to find a stable 
vacuum with a vanishing cosmological constant. On the other hand, it also specifies that 
under special circumstances, such vacua may exist. A particularly illuminating example is 
provided by supersymmetry. 

2.4 Supersymmetry 

Supersymmetry [16] has been used to solve the hierarchy problem in particle physics, 
i.e. it stabilises the mass of the Higgs boson under quantum corrections. Of course, 
with a certain twist of fate, the same type of argument can be applied to show that a 
vanishing cosmological constant is stable under quantum corrections in super symmetric 
models. Supersymmetric models postulate that each particle has a superpartner called 
a sparticle. Fermions of spin 1/2 have scalar partners while gauge fields have fermionic 
partners. The supersymmetry algebra is simply given by the anticommutation relation 

{Q a ,Qa} = 2a^ (2.20) 
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where Q a are fermionic generators and their complex conjugates Qa (they are Weyl spinors 
with a = 1, 2), is the momentum operator and cr M are the three Pauli matrices a' 1 and 
a = 1. The vacuum is represented by a state vector |0 > which is invariant under 
supersymmetry, i.e. it is annihilated by 

g«|0 >= (2.21) 

Specifying \i = and taking the vacuum expectation value of the anticommutation relation 
leads to 

<0|P |0>=0 (2.22) 

where we have used the fact that the vacuum is translation invariant Pj|0 >= 0. Of course 
Po = H is the energy operator (the Hamiltonian) , and therefore we conclude that 

PA = (2.23) 

Strong non-renormalisation theorems exist in supersymmetry which guarantee that this is 
true to all order in perturbation theory. In a nutshell this follows from the existence of as 
many bosons and fermions of equal masses. Loop corrections due to bosons are exactly 
canceled by their fermionic counterparts. 

As a result supersymmetric theories are strongly motivated candidates to explain cos- 
mic acceleration. Let us briefly summarise how they can be built. We focus on the scalar 
sector only. The theories are specified by two functions W{<f> 1 ) where the scalar fields are 
complex and K(<f> 1 , <p % ). The superpotential W and the Kahler potential K can be combined 
to give the scalar potential 

V F = K^FiFj (2.24) 
where is the inverse matrix of Kq = didjK and 

Ft = 3iW (2.25) 

To simplify the discussion, we assume that all the scalar fields are neutral. The kinetic 
terms are simply 

Adn = Kfjd^d^ (2.26) 

Notice that the potential is always positive Vp > and that supersymmetric vacua are 
obtained when F{ = 0. This implies that the cosmological constant vanishes. 

Unfortunately supersymmetry is not a symmetry of nature. Indeed no superpartner 
has ever been observed implying that there must be a mass splitting between superpartners. 
This can be achieved when supersymmetry is broken. Denoting by Msusy the supersym- 
metry breaking scale, the discrepancy between the superpartners leads to a contribution 
to the cosmological constant 

5/>a|susy ~ Mf USY £ 2 (2.27) 

where E is the largest energy scale described by the supersymmetric model. As E > 
-^SUSY > Mz, we find that the breaking of supersymmetry reintroduces a serious fine 
tuning issue. 
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Moreover, supersymmetry must be supplemented with a gravitational sector in order 
to incorporate the Einstein-Hilbert action of general relativity. In this case, the potential 
becomes 

V F = K^FiFj - ?,ml /2 ml x (2.28) 

where 

p. = e 4K/2( d . w + k ^q.k)W) (2.29) 

and the gravitino mass is 

m V2 = e^ K/2 \W\ (2.30) 

The gravitino is the spin 3/2 superpartner of the graviton. In supergravity the potential 
is not positive anymore. Moreover, supersymmetric vacua still satisfy Fj = and have a 
negative cosmological constant 

PA = —3 < m\/ 2 > m P1 (2-31) 

leading to an anti de Sitter geometry. Broken supergravity can have a vanishing cosmologi- 
cal constant provided one fine-tunes the gravitino mass to balance the F-term contribution. 

3. Dark Energy 

We have just seen that symmetry arguments are not enough to guarantee the existence 
of a vacuum with an almost vanishing cosmological constant. It could well be that the 
universe is not in a vacuum state at all and has a dynamical evolution. In this case the 
energy density responsible for the acceleration of the universe need not be constant. This 
evolving energy density has been called dark energy. 

3.1 A fluid approach 

The cosmological constant can be viewed as a fluid pervading the entire universe and 
possessing an equation of state 

w = - (3.1) 
P 

which is constant and equal to -1. In fact, cosmic acceleration does not require the equation 
of state to be —1. We shall first consider fluids for which the equation of state w is constant. 
Accelerating happens as soon as 

w<~ (3.2) 

A natural restriction on the possible value of the equation of states arises from the weak 
energy condition in general relativity: T^ u t^t u > for any time-like vector t^. In other 
words, the energy momentum T^ v must be positive along any particle trajectory described 
by a time- like tangential vector #*. The weak energy condition implies that 

P> 0, w > -1 (3.3) 

where the energy momentum tensor of a perfect fluid is 

T^ v = {p + P)u^u u + pg^ (3.4) 
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and m m is the velocity four- vector of the fluid. Conservation of matter D^T^ V = provides 
the equation 

P = -3H(p + p) (3.5) 
Solving the conservation equation gives the evolution 

P = (-r 3{1+w) P0 (3.6) 
«o 

Using the Friedmann equation for a flat model, one finds 

a = (l) 2 /3(i+-) ao ( 3. 7) 
to 

Of course, this is only valid when dark energy dominates. 

Nothing prevents the equation of state to depend on the energy density too w(p). 
These phenomenological models lead to a time-evolution of the equation of state. The 
simplest of these models is the Chaplygin gas [17] with 

A 2 

P = ~— (3.8) 
The conservation equation can be solved and gives 



a2 + 5 (3 - 9) 

Early in the universe when a is small p ~ B /a 3 like for cold dark matter, while at late time 
p ~ A, i.e. a pure cosmological constant. The emergence of the cosmological constant at 
late time could explain why acceleration is only a very recent phenomenon. Unfortunately, 
this model is under pressure due to its poor properties in explaining structure formation. 
Generalisations with p = —A 2 /p a are acceptable provided a < 0.2 [18]. 

Models with an equation of state w < — 1 violate the weak energy condition and possess 
ghosts. This leads to a catastrophic instability as the scale factor behaves like 

a(t) = (1 - _L) 2 /3(i+™) ( 3 . 10 ) 

^end 

The universe reaches a singularity in a finite time i en d where the Hubble rate diverges: the 
big rip [19]. 

3.2 Scalar field models 

The fluid approach to dark energy is phenomenological. A more fundamental description 
which could be unified with the standard model of particle physics can be obtained using 
a quantum field theory approach. Assuming that Lorentz invariance is not broken in the 
quantum field theory describing dark energy, the fact that the dark energy density has a 
time dependence implies that it must have a spatial dependence too. The simplest models 
for such a space-time dependent energy density are scalar field theories. In the following, 
we shall concentrate on scalar models of dark energy. 
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For simplicity sake, we focus on single field models 



S = J d 4 xV=gl^ - d^<p - V{tj>)\ (3.11) 



The equation of motion for the scalar field in a FRW background is the Klein-Gordon 
equation 

+ 3^ + ^ = (3.12) 

Notice that the expansion of the universe acts as a friction term: the Hubble friction. 
Single scalar field models can be rephrased into a fluid picture where the energy density 
and the pressure are given by 

Notice that the equation of state is always — 1 < w<f, < 1. The dynamics of dark energy 
can be nicely classified in two categories: thawing and freezing models [10]. 

Thawing models are such that the field starts at a non-zero value of the potential and 
stays there due to the Hubble friction. When the Hubble friction is small enough, the field 
starts moving and rolls down towards lower values of the potential. Several possibilities 
can be envisaged. If the potential possesses minima, the field will converge towards one 
of these minima. Depending on the sign of the potential at the minimum where the field 
asymptotically settles down, the universe becomes either a de Sitter space with an eternally 
accelerated evolution if V m i n > 0, or flat space with no cosmological constant if V m [ n = 0. 
If the minimum has negative energy, the universe ends up in a big crunch singularity. 
Phenomenologically potentials leading to V m i D = are favoured. Freezing models are such 
that the potential has no minimum. The field keeps rolling down the potential slope in a 
decelerating fashion until it virtually stops under the effect of the Hubble friction. 

In both cases, the mass of the scalar field needs to be related to the Hubble rate. 
Indeed, if the mass of the scalar field identified with 

< = % < 314 > 

is large compared to H, the field rolls very fast along the potential. The kinetic terms 
dominate and the equation of state is close to w$ ~ 1. On the contrary, if the mass is very 
small compared to the Hubble rate, the field is in a slow roll state. The kinetic energy is 
negligible and uu ~ —1, i.e. the model mimics a pure cosmological constant. The only 
possibility to obtain a dark energy behaviour with an evolution of the energy density is to 
satisfy 

7710 £3 H (3.15) 

This is a strong constraint on the type of dark energy models as it states that the mass of 
the scalar field now must be tiny « Hq ~ 10 -43 GeV. 
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3.3 Attracting and tracking 

In this section we will describe two families of very useful potentials leading to dark energy 
of the freezing type. Thawing models will be exemplified later. The first family consists of 
inverse power law potentials [20] 

M 4+n 

V = — — (3.16) 

The behaviour of dark energy with such a Ratra-Peebles potential is interesting as the long 
time behaviour of the scalar field is largely independent of the initial conditions. Indeed, 
let us consider the radiation or the matter eras where the dominant fluid has an equation 
of state Wb = 1/3 or wb = 0. In this regime, the Klein-Gordon equation has an exact 
solution 

(ftr^t'^ (3.17) 

corresponding to an energy density 

2n 

P4, ~ t ™+ 2 (3.18) 

Moreover this solution is an attractor. Choosing any initial solution for (ft in the early 
universe, the solution always converges to this attracting solution. 
During this regime the energy density decays as 

p<p ~ a~ A 2+™ (3.19) 
This corresponds to an equation of state 

nw B ~ 2 

w* = (3.20) 

Notice that during the matter era, the energy density decays like a - 3n /( n + 2 ) which is slower 
than the matter evolution in a~ 3 . This implies that the dark energy density eventually 
catches up with the matter density and becomes dominant. When this is the case, the 
solution of the Klein-Gordon equation leaves the attractor (only valid when matter or 
radiation dominate) and the universe starts accelerating. One of the nice features of this 
model is the independence of initial conditions. Dark energy becomes always dominant at 
some point. 

The value of the field now is also easily determined noticing that the mass ~ 
V((ft)/(ft 2 . Using the fact that V((fto) ~ p c o and ~ Hq, we find that 

(fto ~ mpi (3-21) 

Hence independently of the initial conditions, the field freezes around a value close to the 
Planck scale now. Of course this can only happen if the scale M in the potential is such 
that 

M 4+n » p c0 m^ (3.22) 

Once the scale M is fixed, the dynamics of the model follow: irrespective of initial conditions 
the field (ft reaches the attractor during both the radiation and matter eras. It leaves the 
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attractor in the recent past and stops at a value close to the Planck scale where acceleration 
starts. 

Another interesting class of models involves exponential potentials [21,22] 

V{(f>) = M 4 exp(-AK 4 0) (3.23) 

where M is an overall scale and A > 0. The potential decreases towards zero at infinity. 
There are two remarkable regimes for these models. First of all if dark energy dominates, 
an attractor solution with 0,^ = 1 can be obtained provided 

\ 2 <3(l + w B ) (3.24) 

In this case, the asymptotic behaviour of the scalar field energy density is such that the 
equation of state is 

W4> = -1 + — (3.25) 

Acceleration occurs when 

A < V2 (3.26) 

In these models, the present state of the universe is transient, matter becomes more and 
more diluted until eventually the energy density of the universe becomes scalar field dom- 
inated. 

For larger values of A 

\ 2 >3(l + w B ) (3.27) 
the scalar field energy density has the same equation of state as the background fluid 

= w B (3.28) 

and the scalar field energy density fraction is fixed 

A, = (3-29) 

i.e. the ratio of the matter energy density over the scalar field energy density is fixed. 
This solution is said to be scaling. The tracking solution is such that the background 
fluid and the scalar field go hand in hand, hence the name tracking. Of course, in this 
regime, acceleration cannot occur. This can be remedied when two exponential terms are 
present [12] 

V{4>) = M 4 (exp(-AiK 4 (<A - <f> c )) + exp(-A 2 K 4 </>)) (3.30) 

Assume that \\ > 3(1 + wb) and A| < 2. As long as <j> < \i4> c /(\\ — A2) the potential 
is dominated by the first exponential. The field converges to the tracking behaviour and 
follows the background evolution. Eventually, the field overcomes 4> > M<t>c/(M — ^2) and 
the second exponential kicks in leading to an accelerated expansion. Of course these models 
suffer from a fine tuning issue as the scale M or equivalently the initial value of 4> must be 
tuned to reach Qj, now. 
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Another class of models is particularly useful. The potential is given by [23] 

V(4>) = V e- XK ^(a + (k 4 - bf) (3.31) 

where A 2 > 3(1+u/b). These Albrecht-Skordis potentials are such that in the early universe 
the exponential behaviour dominates and the scalar field tracks the background fluid. Later 
the field feels the presence of a local minimum and gets trapped there mimicking the role 
of a pure cosmological constant. 

For more general potentials, there is a very useful criterion allowing one to distinguish 
models leading to the asymptotic domination by the scalar field and therefore cosmic 
acceleration [24] . Defining x = -^=^ , the effective coupling constant 

dlnV 

A * = -^sr < 3 - 32 > 

satisfies the differential equation 

^± = -^6\Z(T-l)Hx (3.33) 
Assuming that the field rolls towards large values, and provided 

r > 1 (3.34) 

where 

VV" 

r = -yn ( 3 - 35 ) 

with ' = d/d(fi, the coupling constant converges towards the fixed point A^ = 0. In this 
case, the slope of the potential vanishes, the scalar field energy dominates eventually and 
acceleration occurs. For instance for Ratra-Peebles potentials 

r = 1 + - > 1 (3.36) 

n 

while T diverges at the minimum for Albrecht-Skordis potentials. 
3.4 Coupled dark energy 

So far dark energy has been completely decoupled from the rest of physics. This is not 
always the case, especially in models derived from string theory or extra dimensions. As a 
simplified first step, let us focus on a direct coupling between dark energy and cold dark 
matter [25]. In this case, both the dark matter energy density and the dark energy one are 
not separately conserved but satisfy 

p m + 3Hp m = a^pmcf) (3.37) 

and 

fa + 3H(l + w^p^ = -a^praij) (3.38) 

These equations will be explained when we deal with scalar-tensor theories. We restrict 
our attention to single exponential models with a coupling A. The behaviours obtained in 
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the decoupled case can be generalised. The scalar field reaches an attractor with = 1 
when 

x («2 +12(1+ ))l/2_ 

A < — (3.39) 



In this case the equation of state is 



A 2 

W4> = -1 + — (3.40) 



with acceleration for A < 2. When A is larger 

A > (0 * + 12(1+ ; B))1/ ' 2 -°* (3.41) 

the scalar field has a scaling behaviour where the ratio of the matter and scalar field energy 
densities is constant. It is convenient to define an effective equation of state 

WeS = , (3.42) 

P<j> + PB 

given by 

w c q = — — (3.43) 

A + 

Acceleration can occur when 

a, > (3.44) 
In this case, the energy density fraction of the scalar field is 

o «0(«0 + A) +3(l + w B ) , 

Q + = WTW ( } 

For large values of a^, one can get Q.^ ~ 0.7. Unfortunately explicit models using this 
attractor suffer from phenomenological problems such as a short matter dominated era. 

3.5 Phantom dark energy 

Scalar field models can easily have an equation of state < —1. Consider the action 

S = j d 4 x^[^ + d^d^ - V {<!>)] (3.46) 
where the sign of the kinetic terms has been reversed. The equation of state becomes 

f - vm 

which can satisfy w<p < — 1. We have already seen that this would lead to a cosmological 
big rip singularity. One can even show that the effective field theory with a wrong sign for 
the kinetic terms cannot be valid for energies higher than 

E < 3MeV (3.48) 

coming from the production of a gamma-ray background according to: vacuum — > 20 + 27 
where the vacuum creates particles from nothing. Of course, it seems difficult to believe 
that such a low energy effective theory with ghosts exists at all [26] . 
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3.6 Dark energy difficulties 

Dark energy suffers from two major drawbacks. A very natural problem follows from the 
late appearance of cosmic acceleration. If this is due to a pure cosmological constant then 
the instant when cosmic acceleration sets in only depends on known features of the radiation 
and matter eras. Indeed the decrease of the matter density with time is known and the 
instant when it becomes subdominant to the cosmological constant can be determined. On 
the other hand, once the energy density of dark energy becomes dynamical and evolves in 
time, one may wonder why acceleration starts so late in the history of the universe. This 
can be reformulated as a question: what is so special about the present epoch that the 
matter density and the dark energy density are of the same order of magnitude. This is the 
coincidence problem. In the Ratra-Peebles class of potentials, such a coincidence requires 
to tune the value of the scale M. For exponential potentials, the initial condition of the 
field </>i n i has to be carefully chosen. Of course, in the Ratra-Peebles case one could argue 
that the scale M should be derived from a more fundamental theory therefore leading to 
a solution of the coincidence problem. So far, no solution to the coincidence problem has 
ever been uncovered. 

Another important difficulty springs from the role of quantum corrections [27]. We 
have already seen that the value of the potential can be shifted by quantum corrections. 
As a result, dark energy is a valid concept provided the cosmological constant problem 
has been solved, i.e. the minimal value of the potential vanishes. In thawing models, this 
requires the existence of a minimum with a vanishing potential. In freezing models, the 
asymptotic value of the potential when the field rolls to very large values has to vanish 
too. We have seen that no intrinsically field theoretical explanation has been given to this 
puzzle so far. It might be that the vanishing of the minimal value of the potential has 
an anthropic origin. If this could be ascertained, dark energy with specific potential could 
lead to an explanation of cosmic acceleration due to a tiny vacuum energy density. 

Unfortunately the dark energy mechanism is also jeopardised by the smallness of the 
dark energy mass « Hq. This is such a tiny scale compared to particle physics ones 
that it can be drastically shifted by quantum corrections due to the decoupling of heavy 
particles of mass M 

R 2 M 4 

S < " Wqjr <"») 

where the coupling f3 = for scalar-tensor theories. The correction term is too large as 
soon as M > 10~ 3 eV corresponding to neutrino masses if the coupling (3 = 0(1). Of 
course corrections to the dark energy mass coming from standard model particles are way 
too large unless 

P < 1(T 40 (3.50) 

This would lead to an almost complete decoupling of dark energy from matter. In all other 
cases, stabilising the mass of the dark energy field is a very hard problem. 

3.7 Pseudo-Goldstone bosons 

Goldstone bosons appear when a global symmetry is broken at a scale /. Typically, consider 
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a U(l) global symmetry with a complex scalar field acquiring a vev < <E> >= / at the 
minimum of its potential V(&). Writing 

$ = f e i<t>/V2f ( 3 _ 51 ) 

the boson ^ is a Goldstone boson with normalised kinetic terms. The Goldstone field has 
no potential and is therefore massless. Moreover the model has a residual shift symmetry 
<j) — > <p + c where c is real corresponding to the original U(l) global symmetry. This implies 
that couplings to fermions of the type 



--int 



$<W>i> (3-52) 



are forbidden. A potential for the Goldstone boson can be obtained by a small breaking of 
the U(l) symmetry 

<j) _|_ $ 

-^-breaking = M (3.53) 

where fj, <C /. This leads to the potential 

V (cf > ) = ^ cos -^j (3.54) 

This model is of the thawing type [13]. Initially the field starts away from the minimum 
of the potential and stays there due to the Hubble friction. Later, it converges to the 
minimum of the potential. Cosmic acceleration implies that 

H ~ 10~ 3 GeV (3.55) 

The mass of the pseudo Goldstone boson is then 

,4 



.2 _ V 

it 

Requiring that m$ « Hq leads to 



m ~ 2/2 ( 3 - 56 ) 



/ « mpi (3.57) 

Corrections to the mass of the pseudo Goldstone boson arise from derivative couplings of 
the field <3? to fermions. The lowest order term is 



Ant = —f3~H> (3-58) 



T 3 

This non-renormalisable operator is suppressed by the Planck scale as the model must be 

valid up to this scale to accommodate / ~ mp\. To leading order, this gives an interaction 
term for the pseudo Goldstone mode 

Ant ~ ( 3 - 59 ) 

The mass correction due to this interaction is of order 

n ml M 2 

5m ^T^ (3 ' 60) 
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Of course this is a tiny correction. Hence pseudo Goldstone boson solve the mass problem 
of dark energy. 

An explicit realisation of the Goldstone idea can be obtained using neutrino physics 
[28]. Indeed consider a family of three right handed neutrinos Ni whose Majorana mass 
matrix is due to the interaction with 6 scalar fields <&jj according to 

C N = -X^ijNiNj (3.61) 

The model is invariant under a U(l) 3 symmetry obtained by changing the phases of the 
right handed neutrinos. When the couplings \j = the model is invariant under U(6) 
corresponding to changing the phases of the six scalars. Now when &ij acquires a vev 
< >= fij the six Goldstone bosons separate into three true Goldstone bosons and 
three pseudo-Goldstone bosons due to the Majorana interaction terms. A potential for the 
pseudo-Goldstone is generated at one loop by the Coleman- Weinberg mechanism 

y = ^tr(MMtMMtln(^)) (3.62) 

where Mij = Xijfije 1 ^^^^ with no summation involved and E is the limit of validity of 
the model. The fields <pij are the six Goldstone bosons. This gives a potential to the three 
pseudo Goldstone bosons of the form 

y(^)« / u 4 cos(-^) (3.63) 

where fi 4 ~ M . Now neutrino oscillation experiments are compatible with n ~ 10~ 3 eV 
providing a possible link between neutrino physics and dark energy. 



4. Scalar-Tensor Theories 

Many well-motivated models lead to scalar-tensor theories where a scalar field couples to 
matter on par with the gravitational field. Another type of models, the f(R) theories, are 
also scalar-tensor theories. 

4.1 Jordan vs Einstein 

The action governing the dynamics of the field 4> in a scalar-tensor theory is of the general 
form 

S = j ^xV^I^-^W) 2 -^)}-/^^^,^). (4-1) 

where g is the determinant of the metric g^ u , R is the Ricci scalar and ip$ are various 
matter fields labeled by i. A key ingredient of the model is the conformal coupling of 4> 

(i) 

with matter particles. More precisely, the excitations of each matter field ipm follow the 
geodesies of a metric g^ v which is related to the Einstein-frame metric g^ v by the conformal 
rescaling 

~g^ = A 2 {<P)g^, (4.2) 



- 21 - 



The Klein Gordon equation is modified due to the coupling of the scalar field to matter 

□0 = -a T+^ (4.3) 

where T is the trace of the energy momentum tensor T^ u and the coupling of the scalar 
field to matter is defined by 

dlllA /, AS 

This is equivalent to the usual Klein-Gordon equation with the effective potential 

V eff (<P) = V(<f ) )-A(<i ) )T (4.5) 

Matter is not conserved anymore. Indeed the fact that the scalar field couples to matter 
implies that 

= a^T. (4.6) 

This can be also seen by looking at the Lagrangian involving matter and the scalar field. 
Let us first examine the effect of the scalar field on couplings and masses. The scalar field 
does not couple to gauge field 

W = -^2 / d 4 x^F^F^. (4.7) 

On the other hand, the scalar field couples to both scalars and fermions with a mass 

m(<f>) = A((f))mo (4.8) 

where mo is the bare mass in the Lagrangian. This leads to the possibility of varying 
particle masses. 

In the cosmological context 

T = -p + 3p (4.9) 

implying that the scalar field does not couple to radiation. Moreover the conservation 
equation becomes 

p + ?,H{p + p) = a 4> ^{p-3p) (4.10) 

In the case of cold dark matter, this is the equation we have used in the coupled dark 
energy case where is constant. The effective potential depends on the amount of matter 

V eff (<l>) = V(<f>) + A(<f>)p m (4.11) 

This has spectacular consequences for chameleon fields. 

Matter follows geodesies which deviate from the ones in ordinary gravity 

d 2 x % , dxi dx k dd> , , 

^ + F ^^ +a ^dx- = ° (412) 

in terms of the proper time r and the Christoffel symbols calculated with the metric g^ u . 
In the Newtonian case where 

ds 2 = -(1 + 2<S> N )dt 2 + (1 - 2$ N )dx i dx i (4.13) 
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and assuming that is slowly varying along the particle trajectory, this reduces to 

(fix 1 

_ = -&{& N + a^4>) (4.14) 

This can be interpreted as the motion of a particle in the effective gravitational potential 
which is 

§ = $at + a^K4^ (4.15) 

This is the potential acting on a single particle. The scalar field induces a modification of 
gravity. 

Modifications of gravity are drastically illustrated for models with V = 0, i.e. a 
massless scalar field and p = moS^ for a massive point particle. In this case the static 
Klein-Gordon equation reads 

A<p = Q K 4 mo(5 (3) (4.16) 

whose solution is 
implying that 

¥ = -(l + 2a})^ (4.18) 
leading to a modification of Newton's law depending on the coupling 

G N , ef f = iG Nl (4.19) 

where 7 is one of the Eddington post-Newtonian parameters 

7 = 1 + 2aJ. (4.20) 

These results have been obtained in the Einstein frame where the Einstein-Hilbert term in 
the action is not modified [29]. 

There is another representation of the same theory in the Jordan frame obtained by 
a Weyl rescaling amounting to considering that g^ v is the dynamical field representing 
gravity. Using the transformation law 

R = A 2 '((/)) (R + 6g^D^D v In A(<f>) - Gg^D^ In A(<f>)D v In A(<f>)) (4.21) 

the Jordan frame action reads 



8 = f A^{^ij-^(l-64)(^-^}-/A£ m («,S„„), 

(4.22) 

Notice that in this frame matter is conserved. The Jordan frame is such that matter 
couples minimally to gravity. The modification of gravity appears due to the A 2 ((f)) in 
the Einstein-Hilbert term. In particular, the measured Newton constant in Cavendish 
experiments is 

G N = (l + 2a 2 )A 2 (<p)G N (4.23) 
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The coupling constant is strongly constrained by the Cassini experiment [30] 

a\ < 10~ 5 (4.24) 

where in most cases A((J)q) 1. If the scalar field does not couple to baryons as in coupled 
dark energy, the constraint is relaxed. 

4.2 Violation of the equivalence principle 

Of course, gravitational experiments are not carried out on microscopic particles but on 
macroscopic objects composed of many atoms. The mass of a particular atom can be 
decomposed as [29] 

"atom - MA QCD +a'(N + Z) + 8' (N - Z) + a 3 a QED £ A A QCD , (4.25) 

where A QCD ~ 180 MeV is the QCD scale, iV the number of neutrons and Z the number 
of protons. The quantity M can be written M = (N + Z) + E QCD / Aqcd where E QCD is 
the strong interaction contribution to the nucleus binding energy. The number E\ is given 
by Ea = Z(Z — 1)/(N + Z) 1 / 3 and the quantity a^a QW A. QC - D EA represents the Coulomb 
interaction of the nucleus where a3a QED ~ 0.77 x 10 -3 . Finally the coefficients 5' and a' 
depend on the constituent masses and can be expressed as 

° = \(m u + m d )(b u + b d ) + ^(C n + C p ) + \m e , (4.26) 

^ (mu _ mdKiu _ 6dH ^ (c „_ Cp) _i m ,, 



where m u ~ 5 MeV, m d ~ 10 MeV. The constants appearing in Eqs. flOg ) and (|42p arc 



given by: b u + 6 d ~ 6, b u - b d ~ 0.5, C p q qed ~ 0.63 MeV, C n a QED 0.13 MeV. This 

implies that cr'/A QCD ~ 3.8 x 10~ 2 and 5'/A QCU ~ 4.2 x 10~ 4 . 

The fact that m VL ,m d and m e are dark energy-dependent quantities imply that the 
coefficients a' and 5', and hence m ATOM , are now ^-dependent quantities. We are now in 
a position to estimate the typical gravitational coupling of an atom. From the previous 
considerations, one obtains 

N + Z d ( a 1 \ N-Z d ( 6' \ /A , 

M Qato „ = — ^ (j-) + — M [j-J (4.28) 

where the variations of a' and 5' read 
d 



d 



°' ) 




Aqcd / 


ZJi QCD 


' ) 




Aqcd / 


= 2A 

QCD 



[(feu - b d ) (m u - m d ) + me] 



(4.29) 



They are proportional to the microscopic coupling au. On the other hand, each coupling 
to each type of atom is different due to the different numbers of protons and neutrons. As 
a result, we find that each atom couples as 

a A = (3Aa<f> (4.30) 
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where /3a can be deduced from the previous formulae 

= N + Z {b u + b d )(m u + m d ) +m e Z - N (b u - b d )(m u - m d ) + m e 
1 M 2A QCD M 2A qcd 1 ' j 

Here the violation of the equivalence principle is measured by 

\a A - a B \ , , 

7] AB = 2 ■ (4.32) 

a A + a B 

where a a is the acceleration of the atom A in the gravitational field of a body E, it is given 
by 

riAB = \PA-pB\PEOt% (4.33) 

Now the present limit on the violation of the equivalence principle is rj < 10~ 12 . For a 
typical atom j3 = 0.1 this leads to a strong bound 

a<t> < 1(T 5 (4.34) 

which is much stronger than the Cassini bound. This bound applies to dark energy models 
as the mass of the dark energy field m<f> ~ Hq is almost zero (more precisely the range of the 
fifth force mediated by the scalar field is of the order of the Hubble horizon, i.e. much larger 
than the distances which are gravitationally probed). As a result, one must explain why the 
coupling of dark energy is so weak. For thawing models of the pseudo-Goldstone type, this 
is natural as the direct coupling to matter is forbidden by a shift symmetry. For freezing 
models with runaway potentials, this appears to be another fine-tuning. Fortunately, a 
class of models known as chameleons do not suffer from these problems. 

4.3 Examples 

Let us give non-trivial examples of scalar-tensor models. Consider the low energy theory 
of the Randall-Sundrum model with two branes and matter on the positive tension brane. 
At low energy, the distance between the branes becomes a dynamical field (j) according to 

d= -Zln(tanh^) (4.35) 
V6 

where l~ 2 is the bulk Anti de Sitter curvature. The coupling to matter is 

A(<f>) = cosh (4.36) 
V6 

tanh^l 

= (4-37) 

When the branes are far apart the coupling is small. It is equal to = 1/a/6 for close 
branes. 

Another interesting class of models arises in supergravity when dark energy belongs 
to a hidden sector 

K = Matter + ^D E (0, W = HWcr + W DE {4>) (4.38) 



implying that [31] 
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The only interactions between the dark energy and matter sectors are gravitational. In 
particular the mass of matter fermions becomes 



m fermi on(0) = e K ^ DE ^^ 2 m fcrmion (4.39) 

The prefactor is nothing but the A(4>) factor in these models. Noticing that the dark energy 
field is not canonically normalised £kin = K^d^cpd^cj), this leads to [32] 

= (4.40) 

For canonically normalised fields -Kde = <P4>, this gives a$ = which is of order one for 
runaway models where 4>o ~ mp\. For moduli fields with -Kde = — 3m,p|ln(K4(0 + (f>)) this 
is a<f, = ^p. In both cases, the models are ruled out unless the chameleon mechanism is at 
play. 

In fact, a natural way of guaranteeing the vanishing of along <f> = (ft taken as the 
dark energy direction is to impose that [32] 

Kde(<^) = K be (^-$) (4.41) 

The model has a shift symmetry <\> — > (p + c characteristic of pseudo-Goldstone bosons or 
axions. If dark energy can be derived from a fundamental theory such as string theory, it 
is likely that the dark energy particle will be an axion [33]. 

4.4 Unimodular gravity 

Unimodular gravity [9] is an interesting modification of general relativity where the cos- 
mological constant is not dynamical anymore, therefore solving the cosmological constant 
problem. The gravity Lagrangian becomes 

S = J d 4 x(^+M 4 ) (4.42) 

where det(— g) = 1. The cosmological constant Mq is non-dynamical as it does not appear 
in the equations of motion. In Einsteinian gravity, the cosmological constant receives quan- 
tum field corrections due to the non-trivial nature of vacuum fluctuations. In unimodular 
gravity, the radiative corrections do not lead to the acceleration of the universe despite the 
possible existence of a large cosmological constant. The cosmological constant can be arbi- 
trary and has no consequence on the dynamics of the universe. Of course, this possibility 
does not explain the observed acceleration of the late time universe. All it provides is a 
mechanism for the irrelevance of the quantum fluctuations of the vacuum. 

In unimodular gravity, scale invariance is explicitly broken as det(— g) = 1. The 
breaking of scale invariance results in the effective appearance of a scale M in the theory. 
This scale can be understood as the Lagrange multiplier enforcing the constraint 



S 



J d 4 x[y^ ^ + M 4 - M\^g- - 1))] (4.43) 
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Using the action ( [4.43 ), the equations of motion give det(— g) = 1 and 



R»u - y -fR = -K 2 4 M 4 g^ (4.44) 

Taking the covariant derivative of the latter equation, the Bianchi identity implies that M 
is a constant. This scale is not dynamical and is simply set by the initial conditions of 
the theory. In the Einstein equations it acts as a pure cosmological constant. Of course 
this M term is independent of the vacuum energy Mq . In fact it is a boundary condition 
which has to be set in the early universe. 

Dark energy can be incorporated in unimodular gravity by defining the following scalar- 
tensor theory 

Sum = J d 4 x(/(^)i? - hij^d^d^ - Vtf)) (4.45) 

where the mass dimensions are such that are [(f) 1 ] = 1, [/] = 2, [hij] = 0, [V] = 4, and 
the mass dimensions coincide with the scale dimensions under the rescaling x^ — > X~ 1 x >l . 
Unimodularity can be enforced using a Lagrange multiplier 

S = J tfx^WWR-hij^d^dvp - W))]-M 4 (v^-1)) (4.46) 

As in the case of pure unimodular gravity, one can show that M is a constant. The 
equations of motions enforce det(— g) = 1 and also two sets of equations which coincide 
with the Einstein and the Klein-Gordon equations of the fully covariant Lagrangian 

C = V^l(f(<P l )R - hi^d^d,,? - V(^) - M 4 ] (4.47) 

Hence the dynamics of the theory are determined by this scalar-tensor theory. 

Homogeneity allows us to single out one field $ = (p 1 and define z % = cp 1 /^ 1 . Therefore 
the coupling functions become 

/(^) = $ 2 /V), h^) = h, l3 {z k ), V{P) = $ A V(z l ) (4.48) 

Vacua of the theory are defined by minima of the potential d^iV = implying that V = 
and d z iV = 0. Notice that the vev of $ is not determined at all: it is a Goldstone direction 
preserving the vacuum of the theory. In fact $ is only a pseudo -Goldstone mode as the 
breaking of scale invariance in unimodular gravity lifts the vacuum energy density and 
gives a potential to p. The pseudo- Goldstone mode is a natural candidate for dark energy. 

Let us write the effective action along the Goldstone direction. Defining a constant 
7 such that hij(< z l ^d^d^cfP = Q^d^d^Q we find that the effective action for the 
Goldstone mode is 

CeS = yf=g[{fo<5> 2 R - &y8/*^ " M A ] (4.49) 

where /(< z a >) = /q. So far, we have not included matter fields. They can easily be 
introduced by adding their action to the previous model where the matter fields ij) m couple 
to the metric g^ v . 
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In order to analyse the gravitational and cosmological properties of the model it is 
easier to go to the Einstein frame with 

A 2 {<$>) = (2/qkIcD 2 )- 1 (4.50) 

In the Einstein frame the action reads 

Ses = J rf 4 ^[^ - 6(7 + f )A 2 ^)d^d^ - M 4 A\n (4.51) 

It is useful to normalise 

$ = c Mv / 7n<iVv / 6(7+/o) (4-52) 

leading to the effective action 

v^o \/^o 

S = y>x^[A_Id MW _ M V 4K4 v 7 ^V^^ (4-53) 

We recognise a coupled dark energy model and we can identify 

A = 4c^ (4.54) 

where 

a * = ifl 1 \ (455) 
V 6 (7 + fo) 

The equation of state in this model is given by 

16q2 

1 + w+ = (4.56) 

The Cassini bound implies that a| must be less than 10~ 5 . Therefore, dark energy in this 
model is not very different from a pure cosmological constant. 

5. Chameleons and f(R) Gravity 
5.1 Chameleons 

We have seen that the very small mass scale of dark energy wreaks havoc when coupled 
to matter as it leads to a strong violation of the equivalence principle if = 0(1). This 
conclusion is modified for certain scalar-tensor theories called chameleon theories [34,35]. 
Recall that in the presence of matter the effective potential of scalar-tensor theories is 

V eff (<l>) = V(il>)+A(if>)p m (5.1) 

where p m is the ambient matter density. In the following we shall assume that V(4>) is a 
decreasing function as befitting freezing models where the dark energy field rolls towards 
large values. An interesting phenomenon occurs when A((f>) is an increasing function: the 
effective potential acquires an energy density minimum satisfying 

^k*. = -^(^min)^ (5.2) 
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As p m increases, the position of the minimum moves to the origin. At the minimum, the 
dark energy field has a mass 

2 <?V d 2 A 

This mass depends on the density p m . In a very dense environment, the mass can be large 
enough to shorten the range of the force mediated by the scalar field below the observable 
threshold. Gravity experiments have tested Newton's law down to a scale of about 0.1 mm. 
The range becomes shorter and therefore unobservable provided > 10~ 3 eV [36]. In a 
sparser environment such as the solar system, the mass is generically too small to evade 
strong deviations of planet trajectories. Another effect must be at play called the thin shell 
effect. 

Consider a large spherical object of radius R and density p c . Outside the body, the 
density is much smaller p^. In this environment, the profile of the scalar field is non-trivial. 
Inside the body, the field is constant at a value 4> c corresponding to the minimum of the 
effective potential with a density p c . The field has a rapid variation over a shell of size AR 
below the surface of the body. Far away from the body the field reaches a constant value 
4>oo corresponding to the density p^. The solution outside the body is independent of the 
details of the potential and reads 



<f>(r) = -(T^)(^—f^)— + <Poo (5.4) 



R 3 -R 3 A Me-™^-^ 



47rmpi J \ i? 3 

where M is the mass of the spherical body and R s = R — AR. Over distances shorter than 
the large range of the scalar field force outside the body, this corresponds to a modification 
of Newton's potential by a factor 



W tr)=2ai \-R^)- <5 ' 5) 

For a point-like object where R — > and no thin-shell, we retrieve the 2a^ correction. For 
larger objects, the correction is small when the size of the shell is small 

(5.6) 



R %a ( j > {(j) 00 )mpi^ N {R)' 
This is the case when the thin shell condition is satisfied 

(0oo - (pc) < 6a^(^ 0O )mp / $ A r( J R). (5.7) 

When a thin shell exists for large bodies, the field is called a chameleon. 

Thin shells do not exist for all couplings and potentials. If both are exponential, a 
thin shell cannot be obtained for planets in the solar system. A class of models with nice 
properties is obtained with 

M 4+n 

V{<V) = M 4 + — (5.8) 

They correspond to Ratra-Peebles potentials corrected by a pure cosmological constant. 
The coupling is chosen to be 

Aty) = e a ^ 4(t> (5.9) 
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Cosmic acceleration requires that Mq ~ 10~ 3 eV. The thin shell condition for test ob- 
jects used in laboratory experiments implies that M ~ Mq. This is a numerical coinci- 
dence which is highly intriguing. We will come back to the type of potentials leading to 
chameleons when we discuss f(R) models. 

Cosmologically, chameleons have a very simple behaviour. To simplify the discussion, 
we assume that = 0(1) and almost constant. Then we can derive a lower bound on the 
chameleon mass 

m 2 > 3c^F 2 (5.10) 

where we assume that <p <C mpi in cosmology. This is always the case in known models. 
The fact that the mass at the minimum is always greater than the Hubble rate implies 
that the minimum is a cosmological attractor. After a sufficiently long time, the chameleon 
converges to the minimum. Notice that the minimum evolves as the matter density de- 
creases. Prior to converging to the minimum, the chameleon field can either undershoot 
or overshoot. Undershooting corresponds to an initial chameleon <pi ^min* 

During this 

period the potential is negligible and the chameleon obeys 

£ + 3/ty= (5.11) 
mpi 

During the radiation era, the trace of the energy momentum tensor is almost zero implying 
that the field is stuck at its initial value. It would only start moving during the matter 
era. This is only part of the story as one can show that every time a particle becomes non- 
relativistic the trace of the energy momentum tensor becomes non-zero during a Hubble 
time. As it happens, the field receives a kick and moves rapidly before grinding to a halt 
again. The cumulative effect of these multiple kicks is a variation 

A(j) « a^mpi (5-12) 

Overshooting is the opposite situation where <pi <C <j) c . In this case, the chameleon starts 
very high along the potential and starts tumbling down very fast. During this free fall 
evolution, the field does not feel the potential and eventually stops due to the Hubble 
friction 

0sto P = 4>i + \J 6f^mpi (5.13) 

where is the initial density fraction of the chameleon. After this the chameleon is in an 
undershoot situation. 

A variation of the chameleon of order implies a variation in the particle masses of 
order 

^*-<V^ (5.14) 

ITlf THpi 

Constraints from Big Bang Nucleosynthesis imply that this variation should not be larger 
than 10% between BBN and now. If the chameleon is overshooting during BBN, its vari- 
ation is bounded by y / ^Sl^mpi leading to a bound on the initial fraction < 10~ 2 . 
Generically, the field does not overshoot during BBN. It stops before BBN and then under- 
shoots. If the chameleon is undershooting during BBN, then it receives a dangerous kick 
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due to the electron becoming non-relativistic and leading to a variation of order 0.3mpi. 
This is excluded. Hence the chameleon must have reached the vicinity of the minimum 
by BBN in which case the electron kick is non-operative. This requires an initial value of 
order < m-p\. This large value of the field can be obtained if the chameleon is initially 
overshooting with ^ < 1/6. If the field is around the minimum during BBN then the 
subsequent evolution of the chameleon is negligible compared to the Planck scale leading 
to no problems with BBN. In conclusion, chameleon fields are generically required to be 
close to the minimum during BBN. 

5.2 f(R) models 

A class of modified gravity models turns out to be equivalent to certain scalar-tensor 
theories [37]. Let us consider the theory defined by 

:f(R) 



I 



d X\f— g- 



2k\ 



(5.15) 



where / is an arbitrary function. These models extend the usual Einstein-Hilbert action 
involving high order terms with more than 2 derivatives. It is more convenient to introduce 
a field A and write ( 

(5.16) 



s = h I dixV ^~ 9 [/(A) + f ' {x){R " A)] 



As long as /"(A) 7^ 0, this leads to R = A and the appropriate f{R) action. The field 
A has no kinetic terms. The action is a particular form of a scalar-tensor theory with 
A 2 = 1/ f'(X). In the Einstein frame, the model becomes 



/ 



j4 / 

d xy/-g 



v(4>) 



where 



and 



/'(A) 



v{4>) 



A/'(A)-/(A) 
2*2(/'(A))* 



The coupling is fixed 



1 

76' 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



As a result, the only way f(R) theories can pass the gravitational tests is to impose that 
they must be chameleon theories in disguise. This gives strong constraints on the shape 
of the potential as it must lead to the thin shell behaviour. Equivalently this restricts the 
possible functions f(R)- 

Let us analyse the possible deviations of f(R) models compared to a pure cosmological 
constant. First of all, there is a very mild constraint coming from the requirement that 
there should exists a matter era which is then followed by an accelerated epoch. This is 
realised provided [38] 

RfRR 



m 



Ir 



(5.21) 
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is less than 0.1. Another set of much more stringent constraints follows from the thin shell 
condition [39]. Defining 

* = /'(A) (5.22) 
The equations of motions derived from Einstein equations are 

_ *H _ = ^(T + 2&V) (5.23) 

d CI o 



and 



A = ^ + 4*V-&^ (5.24) 
a 4 a 6 <I> 

interpreted as the Klein-Gordon and the Friedmann equations in conformal time 

ds 2 = a 2 {r])(-dr] 2 + dx 2 ) (5.25) 

We will derive a constraint on the apparent equation of state from these equations. 

When observations are carried out, it is generally assumed that particle physics pro- 
cesses which happened in far away objects are identical to the ones we observe on earth. 
This selects the Jordan frame as the physical frame to interpret data. The equation of 
state is then inferred by identification in the Friedmann equation assuming that the dark 
energy fluid is conserved. This amounts to writing 

a 2 

h2 = ^ = 4^o(p m + Pbe) (5.26) 

where 

4 cS = £ (5-27) 

is the apparent gravitational coupling. This apparent Friedmann equation defines pde- 
Conservation of matter defines the dark energy equation of state 

Pde,?? = ~ 3 "^( 1 + w de)pde (5.28) 
With these identifications the equation of state can be obtained 

(1 + WE)fiDE= _^- 2 ^_ + ( |- 1)!ira (5.29) 

The right hand side is of order of the variation of <j) in Hubble time units. Denoting by A(p 
the time variation of (/) during one Hubble time, we find 

|(1 + w;de)^de| =0(^^) (5.30) 
mpi 

Imposing that objects such as the sun have a thin shell leads to an upper bound on the 
spatial variation of the scalar field 

< % (5-31) 



7Bpi 3 
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where U is Newton's potential of the inhomogeneity to which the object belongs (typically 
a galaxy cluster). Using the fact that the field at infinity increases with time 4> c < (f>(t) < 
4>(to), we find that A</> = <^oo(*o) — 4>oo{t) < 4>oo{to) — 4> c and therefore 

|(1 + u; D e)^de| < | (5.32) 

The largest structures in the universe have U = 0(1CU 4 ) leading to an upper bound 

|(1 + w;de)^de| < 10~ 4 (5.33) 

This is a loose bound which can be supplemented by explicit calculations using constraints 
from laboratory gravitational tests such as the Eotwash experiment. In all known cases, 
the equation of state is bounded and lies extremely close to —1. The fact that f(R) theories 
cannot be distinguished from a pure cosmological constant at the background level is nicely 
counterbalanced by the fact that perturbation theory and the growth of structures are very 
different in chameleon theories (hence for viable f(R) models) 

5.3 Structure formation 

It is actually easy to analyse the growth of structures in f(R) models using the equivalent 
scalar-tensor theory. In this section we will present results which are valid for all chameleon 
models. We assume that the background is such that the chameleon sits at the minimum 
of its effective potential. As we have seen, this must be the case since BBN and therefore 
is the case of interest when structures form in the matter era. We also assume that the 
chameleon energy density is subdominant, i.e. we take z > 0(1). In the linear regime, we 
are interested in the density contrast 




(5.34) 



in the Jordan frame. It is convenient to analyse the perturbation equations in the Einstein 
frame where p = A(<p)p implying that 

5 = 5 — K^a^54> (5.35) 

where 5 = In the following we will calculate 5 and 5<p. 

In the Einstein frame, the absence of anisotropic stress implies that the metric can be 
written 

ds 2 E = a 2 (-(l + 2<S> N )drj 2 + (1 - 2$> N )dx 2 ) (5.36) 
in conformal time. Using ds 2 j = A 2 (4>)ds 2 E and writing 

ds 2 j = a 2 ((-(l + 2^)d7] 2 + (1 - 2$)dx 2 ) (5.37) 

we find that 

§ = ^ N + ^S(/), $ = <!> N --^5(t> (5.38) 
mpi topi 
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The Klein-Gordon equation becomes a relation between the density contrast and the field 
fluctuation 



' k 2 + a 2 m 2 . 



in momentum space. We deduce that 



(5.39) 



(5.40) 



In general, the correction term is small. 
Using the Poisson equation 



we find that 



and 



N 



3a 2 H 2 
2k 2 ' 



a 2 m z v <t>> 



1 + 

1 n^-m z 



$ = - 



1 + 



1 + 4^ 

1 n^<rri z 



3H 2 a 2 5 
2k 2 



3H 2 a 2 5 
2k 2 



This leads to the slip function [40] 



(5.41) 



(5.42) 



(5.43) 



7 ^ 



1 + ^(1 + 2a 2 



(5.44) 



Similarly one can identify a deviation from Poisson's law 

3a 2 H 2 



-/3t<5 



(5.45) 



with 



i + -^(i + 2«2: 



~9 7 



1 + 4^ 



1 



(1 + a 



2 3a 2 g 2 
1 k 2 +d- 



2^T 



(5.46) 



It is clear from these expressions that in chameleon theories, the range of the chameleon 
force 



777,0 



(5.47) 



is crucial. On large scales k/a <C m^,, the deviations from general relativity disappear 
7 = 1, {1=1. This is not the case on scales smaller than A c , i.e. k/a^ where 



1 " 2a l 9 
7 = -— 4, A = l + 2a 



1 + 2a 



(5.48) 
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Notice that \x is tantamount to the fact that Newton's constant for small objects is larger. 
For f(R) theories, we find that 

14 , 
7=5' ^=3 (5-49) 

These are large deviations on small scales. 

In fact, these modifications can be summarised in the equation governing 5. First 
define the divergence of the velocity fluid 6 = diV 1 satisfying the Euler equation 

e' + ne = k 2 ip (5.50) 

where the conservation equation gives 

5' = -9 (5.51) 

Using these two equations we find: 

5"+H6'-l-H 2 ^^6 = (5.52) 
in conformal time and in the linear regime. The scale dependent Newton constant is 

G N (k) = (l+ 2a J m2 )G N (5.53) 

It interpolates between Gn on large scales k/a <C and GN.cff = (1 + 2c^)Gat on small 
scales k/a 3> m^. When the coupling is constant, this equation can be integrated. 
When k/a <C m^, the growth is similar to the one in general relativity 

5 ~ a (5.54) 

whereas gravity is modified on small scales leading to an anomalous growth 

5 ~ a% (5.55) 

where 

-1+ /l + 24(1 + 2^) 

v = V _ _ (5.56) 

This implies that structures grow faster on small scales. For f(R) models, we find that 

-1 + ^33 

v = (5.57) 

Using this result, one can evaluate the time dependence of ^ and <!>. On top of being 
different, Newton's potentials ^ and $ are now time dependent. This is important as 
this leads to the existence of an integrated Sachs- Wolfe effect in the Cosmic Microwave 
Background spectrum. The fact that peculiar velocities are sensitive to ^ while weak 
lensing is sensitive to ^ + $ may also be used to disentangle the possible existence of a 
chameleon description of dark energy (of the f(R) type for instance). Another important 
consequence is the fact that gravitational lensing is sensitive to ^ + hence leading to an 
effect on weak lensing. 
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6. Modified Gravity 



It is plausible that gravity may be modified on large scales where it has never been tested 
experimentally. Such a modification could lead to the acceleration of the universe. 

6.1 Massive Gravity 

The first and most natural modification of Einstein gravity is massive gravity. This is the 
gravitational analogue of Fermi's theory of massive gauge bosons. Like for spin 1 fields, 
the mass term breaks gauge invariance which is general covariance here 

9ixv -> + + d u £n (6.1) 
It is simpler to study massive gravity in a Minkowski background by expanding 

9(iu = ri^u + /V ( 6 - 2 ) 
where h^ v is the graviton field. The Lagrangian of massive gravity reads [41] 
77? 2 r r 

s = -^ U A x [(d x h^d x h^ - 2d IM h^d x h x u + 2d^d u h\ - d^d^) + o^r + b{h»f 

(6.3) 

The first term is the expansion of the Einstein-Hilbert term up to second order while a and 
b parameterise the possible mass terms. Notice that we have explicitly respected Lorentz 
invariance. 

The analysis of ghosts in this model is easier using the Stuckelberg formalism. Define 

= h^v + d^u + %4>ii (6-4) 

where under gauge variations 

V/ -> + + d v ^, ->• - ^ (6.5) 

The vector (j)^ is the analogue of the longitudinal polarisation for massive vector fields. 
This trick restores gauge invariance even in massive gravity. Now the Einstein-Hilbert part 
of the action is only a function of h^ v while the mass terms lead to kinetic terms for (p^. 
Considering only these kinetic terms, we find 

o 

) > I 

^Stuckelberg = [ad^uW + (a + 2b){d^f] (6.6) 

This Lagrangian has no ghost provided it can be written up like in electrodynamics involv- 
ing only F^ u = d^(j) u — d v (j)^. This is possible if 

a = -b (6.7) 

in which case 

o 

^Stuckelberg = —^F^F"" (6.8) 
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where we have introduce the mass 

mg = bnipi (6.9) 

which corresponds to the mass of the transverse traceless part of the graviton. In summary, 
massive gravity has no ghost provided provided the mass term is of the special form 

2 

^massive = ^ (V^ " (K) 2 ) ( 6 - 10 ) 

called the Fierz-Pauli mass term. 

A massless graviton has two polarisations while a massive graviton has five polari- 
sations. The gravitational field emitted by source of energy momentum tensor T^ u is in 
general relativity 

1 



while in massive gravity 



_ 87rGjy 



8ttG n 



V- v ~ 9 , 2 



(6.11) 



(6.12) 



Notice that the main difference is not the pole structure but the appearance of a 1/3 instead 
of a 1/2. This implies that even in the vanishing limit — > 0, massive gravity is not 
equivalent to massless gravity. This is due to the extra polarisations which are present in 
massive gravity. This major flaw is called the van Dam- Veltman- Zakharov discontinuity. 
At the non-linear level for a body of mass M, it has been argued by Vainshtein that there 
is radius 

*v=(-3— H (6.13) 

below which massive gravity becomes non-linear and the vDVZ discontinuity is not appli- 
cable. Unfortunately, the ghost which was suppressed in a Minkowski background using the 
Fierz-Pauli choice does not decouple anymore in a curved background, hence in cosmology. 

6.2 DGP gravity 

Dvali-Gabadadze-Porrati gravity [8] is a modification of gravity at long distance using a 
5d setting. It bears close resemblance with massive gravity as we will see. The 5d action is 



S^^ + S d ^h (6 - 14) 



2nj 

where g a b is the 5d metric and gffi the induced metric on a 4d brane where matter lives. 
This model is particularly popular as it leads to a naturally accelerating brane cosmology. 
Consider the dynamics of the induced metric on the brane which we take to be of the flat 
FRW type. Using techniques which will not be described here, the Friedmann equation is 
modified to [42] 

H 2 ± — = (6.15) 
r 3 
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where matter lives on the brane and 8ttGn = n\. The typical scale 



r = ^ (6.16) 

is the new ingredient of the model. The branch of solution with a + in the Friedmann 
equation is such that, for a low density p <C we retrieve the usual Friedmann 

equation. On the contrary when the sign is — called the self-accelerating branch, the 
Hubble rate tends to 

H -»• — (6.17) 

ro 

asymptotically. Hence the cosmological constant leading to the acceleration of the universe 
is determined by the ratio of the 5d over the 4d gravitational couplings. 

This success is marred by the properties of DGP gravity in the presence of a localised 
source on the brane 



87tGjv 

h 



1 ' rn 



(6.18) 



At small scale, gravity is 4d and the interaction potential is in 1/r while at large scale 
r > ro gravity is 5d with a potential in 1/r 2 . Hence gravity is modified at large scale by 
the opening of an extra dimension. Unfortunately, the 1/3 term reminds us that gravity 
in the DGP model has five polarisations as a massless graviton in 5d. We have seen that 
ro is of the order of the Hubble scale now in order to accommodate the acceleration of the 
universe. Hence gravity would have a wrong tensorial structure in the solar system and the 
DGP model would be ruled out. Fortunately, like in massive gravity, there is a Vainshtein 
scale 

r* = (r S r 2 ) 1/3 , (6.19) 

where r$ = 2GnM is the Schwarzchild radius of the source, below which DGP gravity 
behaves similarly to general relativity. On the other hand, DGP gravity suffers from 
the presence of ghosts on the self-accelerating branch preventing its use as a possible 
explanation of cosmic acceleration [43]. However, the DGP modification of gravity can 
be embedded in string theory [44] with six extra dimensions. This provides an explicit 
regularisation of DGP gravity and cures it from the presence of ghosts. 



6.3 Ostrogradski's theorem 

We have just described two models of modified gravity suffering from the presence of ghosts. 
Gravity could be further modified by adding an arbitrary number of terms to the Einstein- 
Hilbert action involving higher order derivatives. In fact, there is a very strong result due 
to Ostrogradski which prevents one from doing so [45]. 

Consider a Lagrangian depending on a simple variable q and let us focus on time 
dependent configurations only. The Lagrangian becomes 



(6.20) 



where is the jth time derivative. The Euler-Lagrange equation is generalised to 

M"*J 8& = ° (6 - 21) 

Assuming that the Lagrangian is non-degenerate, one can express the Euler-Lagrange equa- 
tion as 

=G(q,...,q( 2N ^) (6.22) 

Solving this equation requires 2N initial conditions. This implies that one must be able to 
define a (2N) dimensional phase space. The purpose of Ostrogradski's construction is to 
exhibit a conserved Hamiltonian in such a phase space, and to show that it is not bounded 
from below. 

Choosing the canonical variables 

one defines the Hamiltonian 

N 

H = J2PiQ (i) -L (6.24) 

i=i 

The Hamiltonian is conserved and generates the time evolution 

r) J-f FT 

* = * = -5a (6 - 25) 

If the system is non-degenerate, one can solve the equation 

dL 



N 



(6.26) 



and express 

q (N) =F(Q 1 ,...,Q N ,P N ). (6.27) 

As a result, the phase space contains 2A^ coordinates (Qi, ■ ■ ■ ,Qn) and (Pi, . . . , P/v). It is 
the clear that the Hamiltonian 

H = P x Qi + ■■■ + Pn-iQn + P N F - L(Q U ...,Q N ,F) (6.28) 

is unbounded from below in the (N— 1) direction (Pi, . . . , P/v-i). This is due to the presence 
of ghosts and is lethal for the theory. A possible way out is to consider degenerate theories 
or gauge theories like ordinary gravity where constraints remove some of the ghost-like 
directions. 

7. Violation Of the Cosmological Principle 

We have just seen that constructing a modified theory of gravity is problematic. As it is 
also very hard to understand the properties of dark energy, one may then doubt the validity 
of the cosmological principle. 
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7.1 Inhomogeneous universe 

The apparent acceleration of the universe could be due to a local inhomogeneity whose 
effect would mimic dark energy [46]. This can be rendered rigorous by using the local 
Einstein equations in a patch around us and averaging over the volume of this patch. The 
averaged evolution of the universe does not only respond to the average density of matter 
as the Einstein equations are non-linear [47] . 

Let us assume that the universe is filled with a pressureless fluid and that it can be 
foliated according to a global time coordinate 

ds 2 = -dt 2 + h ab dx a dx b (7.1) 

where x a are spatial coordinates. The expansion rate of the universe becomes a tensor 

e a b = \h ac k ch (7.2) 

Let us define = and the shear tensor = 6£ — The Einstein equations give 
the Friedmann equation 

R {3) + \® 2 - 2a 2 = 16irG NP (7.3) 
where a 2 = a a !,a ab and the equation 

#(3) + @2 + @ = i 27T g nP (7.4) 

where is the curvature of h ab . These equations can be combined to give the Raychaud- 
hury equation 

e + -e 2 + 2(j 2 + 47rGjvp = (7.5) 

to describe the local geometry of the patch. 

To describe global quantities, one needs to take the spatial average over a comoving 
volume D 

< A> D = — I d 3 xVhA (7.6) 
Vd Jd 

where Vd = J D d 3 xVh. It is convenient to define the global scale factor 

al = ^- (7.7) 
compared to some initial volume. Using the commutation relation, 

d< ^ t >D -<A > D =< AQ > D - < A > D < 9 > D (7.8) 
the averaged Friedmann and Raychaudhury equations become 

Hl = ^ L <P>D ~\{Qd+ < i? (3) >d) (7.9) 

and -An 1 

a D 47tGat 1 , 7in , 

— = 5— <P>d+i:Qd (7.10) 

ao 6 6 
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where Hp = ap/ciD and 



Qd = \(<Q 2 >d- <Q> 2 D )-2a 2 . (7.11) 



As expected, the evolution equations for the averaged cosmology receive corrections terms 
due to the non-linearities of the Einstein equations. Surprisingly, these equations are 
equivalent to the ones obtained with a dark energy fluid 4>D whose potential reads [48] 



and 



12 tti± \ Qd 



° = v ' (fe) - ^ (7 ' 13) 

Of course this result is very intriguing as it states that dark energy could be the result of 
the local dynamics of our local patch of universe. To go further, one needs to calculate the 
averages in a particular model. 

7.2 Tolman-Bondi universe 

The CMB gives a very precise indication that the universe is isotropic at large. If accel- 
eration is due to a local inhomogeneity in our surrounding patch, the phenomenological 
constraints coming from SN la and the CMB are quite tight [49] . First of all, outside our 
patch, the Hubble rate should be as low as h ~ 0.45. Moreover one must accomodate a 
patch reaching as far as the first supernovae z ~ 0.05 implying the existence of a large void 
of size between 160//iMpc and 250 /h Mpc. Inside this underdensity, the density contrast 
must be as large as 5 ~ —0.4. All these constraints can be met using special solutions of 
Einstein equations based on the metric 

ds 2 = -dt 2 + ^^/iK dr 2 + R 2 (r, t)(dd 2 + sin 2 6dcp 2 ) (7.14) 
1 — k(r)r 2 

known as a Tolman-Bondi universe. Notice that these space-times are isotropic but inho- 
mogeneous. They generalise the FRW solutions and represent an explicit violation of the 
Copernican principle: the observers have a special place in the universe. The Friedmann 
equation can be expressed in terms of 

~ Hr)r 2 = R(r,tf (7.15) 

where the comoving mass function M(r) inside a comoving patch of radius r is time inde- 
pendent and is related to the density inside the patch 

M(r) = 471-/" dr'p(r',t)R'(r',t)R 2 (r',t) (7.16) 
Jo 

The function k(r) is the local curvature inside the patch. Let us consider a patch where 
the curvature is negative. The solution of the Friedmann equation is 

R(r,t) = - G " Mir 2 \ cosh v -l) (7.17) 
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where time is given by 



t ~ to(r) = GnM( "1 (sinh V - n) (7.18) 
(— k(r)r z ) ' 

The Tolman-Bondi universes have two arbitrary functions: the curvature k(r) and the 
initial time to( r ) interpreted as the time of the big bang at the location r. The volume 
inside the comoving sphere is given by 

r B = fc r«» ail.) 

jo a/ 1 — k{r)r z 

It has been shown that the averaged acceleration parameter [50] 

d D a D ,_ on x 

to = ~r~2 (7.20) 

is negative when 

hr- > ^ (7-21) 

where 

jo VI — kr z Jo yl — Ar z 

and 

I kr 2 = 4tt f D dr^£EL (7.23) 
Jo V 1 — A;r 2 

In particular choosing 

k(r) = - T ^, < a < 2 (7.24) 

leads to an effective acceleration. It would be very interesting to have more general results 
and a correspondence with the equivalent dark energy description. 



8. Conclusions 



These lecture notes have presented some of the issues related to the observation of cosmic 
acceleration. The list of topics which have been covered is rather small compared to the 
breadth of the existing literature. Four points of view have been commonly advocated. 
The first one which is favoured by some astronomers is that a pure cosmological constant 
is enough to explain all the physics behind cosmic acceleration. I hope I have managed to 
convey the idea that a sheer cosmological constant is not entirely satisfactory and calls for 
a deeper explanation. One of the possibilities could be that the value of the cosmological 
constant has an anthropic origin. Another explanation could be dark energy although 
it requires a huge cancelation too and explains only the observation of a tiny residual 
cosmological constant. In the class of dark energy models, axion-like models are the best 
motivated ones from the particle physics point of view. Another alternative could be 
modified gravity. As I have recalled, this is a tortuous route strewn with difficulties such 
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as the existence of ghosts. In view of its relevance, the challenge of building an appropriate 
infrared modification of general relativity is worth pursuing. Finally the violation of the 
Copernican principle and the role of inhomogeneities are very hot topics which may hold 
the secret of cosmic acceleration. 
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